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The dynamics of a system composed of two coupled optical cavities, each containing a single 
two-level atom, is studied over a wide range of detuning and coupling values. A description of the 
field in terms of delocalized modes reveals that the detuning between the atoms and these modes 
is controlled by the coupling between the cavities; this detuning in turn governs the nature of the 
dynamics. If the atoms are highly detuned from both delocalized field modes, the dynamics becomes 
dispersive and an excitation may be transferred from the first atom to the second without populating 
the field. In the case of resonance between the atoms and one of the delocalized modes, state transfer 
between the atoms requires intermediate excitation of the field. Thus the interaction between the 
' two atoms can be controlled by adjusting the coupling between the cavities. 
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I. INTRODUCTION 



Cavity quantum electrodynamics (QED) provides a natural setting for distributed quantum information processing 
(QIP) [1]. Atoms in highly detuned cavities may be used to store qubits in decoherence-free environments, while 
photons are well suited to distribute information throughout the system. The strong coupling that can be achieved 
I— I between atoms and cavity field modes can be used to perform logic operations 0, Q • The physics of cavity QED 
i-C ' can also be realized in a solid-state analog known as "circuit QED," in which superconducting qubits are coupled to 
^p-i; stripline resonators @|- Circuit QED offers several advantages over traditional cavity QED, including extremely 
high dipole moments that result in large couplings, a fixed number of qubits in a cavity, and improved scalability. 
One requirement of distributed QIP is the coupling of distant qubits in order to perform state transfer, entanglement 
' generation, or quantum gate operations between separate nodes of the system. In cavity QED, two atoms may be 
coupled by connecting the cavities in which they are confined via an optical fiber. In an early proposal 0] for state 
transfer in this system, the atomic state is converted into a superposition of Fock states of the cavity field mode. 
This photonic state then propagates along the fiber to the second cavity, where it is transferred to the second atom. 
Such a technique, however, is susceptible to decoherence due to losses from the cavity or fiber modes. To decrease 
' these losses, systems of adiabatic transfer have been devised, using a sequence of laser pulses to transfer an atomic 
00 '. state without populating intermediate states. It is possible to minimize the population of the cavities 3 O'" the 
OO ■ fiber . Alternatively, the fiber mode can be eliminated from the dynamics if the cavity-fiber coupling is sufRciently 
strong [lO| . High-fidelity swap and entangling gates may then be produced by turning the atom-cavity interaction on 
■r^ljj- i and off. All these proposals, however, still involve excitation of the field in either the cavities or the fiber. 

■ Field excitation can be avoided entirely by coupling two qubits to a single field mode. In cavity QED, this requires 
OO , two atoms to be confined within the same cavity illi] : in circuit QED, two superconducting qubits are coupled to one 
' stripline resonator If both qubits are highly detuned from the field mode, the interaction with the field mode 
becomes dispersive. Thus the field mode is effectively eliminated from the qubit dynamics and an excitation may 
pass between the two qubits without ever populating the field mode. This interaction can produce decoherence-free 
. , state transfer and entanglement generation. However, some practical difficulties remain with this scheme. First, 
i the two qubits, interacting with the same field mode, may not be sufficiently independent to act as separate nodes 
" " ■ in distributed QIP. Second, there is no way to regulate the qubit-qubit coupling without manipulating the qubits 
themselves. 

Recently, considerable theoretical effort has been devoted to a class of coupled-cavity models that promise to 
overcome the pro blem of individual adressability. Inspired by experimental advances in photonic crystals [IjI , optical 
microcavities |l4j, and superconducting devices [3], these models typically describe a series of optical cavities, each 
containing one or more atoms or qubits. Photons are permitted to hop between the cavities. Much of the work 
on coupled-cavity models has focused on similarities to the Bose-Hubbard mode l [l5| and the prospect of observing 
quantum phase transitions between Mott insulator and superfluid states p^. [TtI [iSl. Il9l. [20l. [2ll. l22l l23j. 

Coupled-cavity models have potential applications in QIP as well, since it should be possible to control and measure 
individual lattice sites. A proposal for using a coupled-cavity system to generate high-dimensional entangled states of 
the field was put forward in [23| . Angelakis et al. [25| and Hartmann et al. [1^ have discussed the creation of cluster 
states for one-way quantum computation. State transfer along an array of polaritonic qubits in a coupled-cavity 
system was studied in [27]. A protocol for producing two-qubit quantum phase gates in a system of quantum dots 
coupled to photonic crystal nanocavities was outlined in [28|. 

Some work has also been done on the specific problem of transferring atomic states between two sites in a coupled- 
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cavity system. Skarja et al. [HI studied a system of two weakly coupled cavities. When atoms are passed successively 
through each cavity, quantum interference effects appear in the final state of the second atom due to the existence 
of two field modes. The authors focused on the case of exact resonance between the atoms and cavity fields. More 
recently, Nohama and Roversi [SO*] addressed the problem of atomic state transfer in a system of two cavities, each 
containing a trapped two-level atom. They showed that perfect state transfer occurs periodically, but the analysis 
was restricted to the case A = —A, where A is the detuning between the atom and cavity field and A is the photon 
hopping strength. 

In the present paper, we examine the dynamics of a simple two-site coupled-cavity model over a large range of 
values of the qubit-cavity detuning and the photon hopping strength, focusing particularly on atomic state transfer. 
Although the exact dynamics of the system is rather complicated in the general case, there are several limiting regimes 
in which it may be simplified. In each of these limits we derive an effective Hamiltonian that captures the essential 
physical processes and study the resulting dynamics. The limits of large hopping and large detuning each lead to 
dispersive interactions, allowing state transfer between the atoms without excitation of the field. When the hopping 
and detuning are approximately equal, a resonant interaction between the atoms and a delocalized mode of the field 
occurs. Fast state transfer involving excitation of the field is possible in this case. We find, then, that the hopping 
parameter provides a flexible source of control over the dynamics of the system. 



II. COUPLED CAVITY SYSTEM 



The system under consideration consists of two coupled cavities, each modeled as a single harmonic oscillator mode 
of frequency w/; each cavity contains a single two- level atom with transition frequency = (jJ/ -I- A [36l |. The atom- 
cavity interaction is described by the Jaynes-Cummings model [3l|, using the rotating-wave approximation ^]. The 
two cavities are coupled so that photons may hop between them. The Hamiltonian governing the system is {h = 1) 
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H ^'^\ijJa\e)j{e\+ujfa^^aj+g{aja.j++a]aj-^ + A {a\a2 + a\a^ , (1) 

where g is the atom-cavity coupling strength and A is the cavity-cavity hopping strength. The operator aj (aj) is the 
lowering (raising) operator for the field in cavity j. The states \g)j and \e)j represent the ground and excited states, 
respectively, of the atom in cavity j. The operator dj- {d'j+) is the atomic lowering (raising) operator for cavity j. 

The first two terms under the summation in Eq. ([T]) correspond to the energy of the bare atoms and field modes, 
respectively. The third term represents the transfer of excitations between the atom and field in each cavity. The 
exchange of photons between cavities is expressed in the final term of Eq. ([1]) . 

It will prove useful to rewrite the Hamiltonian in terms of delocalized field and atomic modes, which are given by 
the antisymmetric and symmetric superpositions of the states at the localized sites. The raising operators of the new 
field modes are given by 



rh\ = -^ (a\ + an ; (2) 



V2 

similarly, the raising operators of the new atomic modes are given by 

Sl+ = - CT2+) , 

S2+ = ((71+ + (72+) . (3) 

In terms of these operators, Eq. ([1]) becomes 

[^tJaSj+Sj- + -I- g (nijSj+ + m]sjj^ . (4) 
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In this form, the delocalized atomic mode sj and the delocalized field mode rrij are governed by a Jaynes-Cummings- 
type interaction. No coupling occurs between the odd and even modes, so that the system Hamiltonian takes the 
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form of a sum over two independent Jaynes-Cummings Hamiltonians. The transition frequency of the delocahzed 
atomic modes remains unchanged at oja- However, the frequencies of the delocahzed field modes are shifted from the 
bare cavity frequency: the frequency of mode rrij becomes ujj , where wi = Uf — A and lj2 — lo f + A. The detuning 
between the atomic mode Sj and the field mode rrij is thus given by Aj , where 



Ai = A 
A2 - A- 



A, 
A. 



(5) 



Our study of the dynamics of the system is focused on the problem of state transfer between the atoms. Atom 1 is 
prepared in an arbitrary superposition of the ground and excited states, while atom 2 is in the ground state and both 
field modes are in the vacuum state: 



|vI/(0)) = (cos9\g), 



-e*^sin0|e)i) |g)2|00). 



(6) 



State transfer between the two atoms occurs when atom 2 is in the same superposition of ground and excited states 
that atom 1 had initially. That is, at some time t, the state of the system is given by 



\^{t)) - (cos^|5)2 + e'*sin(?|e)2) |00). 



(7) 



Since the Hamiltonian conserves total excitation number, the dynamics is restricted to the subspace containing zero or 
one excitations. We can thus write the state of the system at time t as — |V'o(^)) + \ ^pl{t)). The zero-excitation 

component \ipQ(t)) ~ \g)i\g)2\00) is invariant under the action of the Hamiltonian, so we need only consider the 
dynamics of the single-excitation component lipiit)). 

A general state (not necessarily normalized to 1) in the single-excitation subspace may be written as 



l^i(t)) =a(i)|.g)i|.g)2|10)+6(t)|5)i|.g)2|01)+c(t)|e)i|5)2|00)+d(t)|5)i|e)2|00). 



(8) 



The coeflScients of \ipi{0)) that correspond to the initial state in Eq. ^ are a(0) = 6(0) = d{0) = and c(0) — e"^ sin 9. 
The conditions for perfect state transfer at time t are, therefore, a(t) = b{t) = c{t) = and d{t) — e**^ sm9. 

The time evolution of the general one-excitation state given by Eq. ([5]) may be calculated exactly from the 
Schrodinger equation. In terms of the delocahzed modes aj{t) = a{t) + {—iyb{t) and (3j{t) = c{t) + {—iyd{t) 
where j = 1, 2, the solutions are given by 



aj{t) 



i/2)t 



a, (0) cos {,y,t) + - ( ^a, (0) ~ g(3, (0) ) sin (i.,t) 
/3,(0) cos (i^.t) - - ( ^/3,(0) + ga,{0) ] sin {,y,t) 



(9) 



where 




(10) 



These equations describe complicated dynamics. However, there are several limiting regimes in which the Hamiltonian 
and the resulting dynamics can be substantially simplified. This is most readily seen in the delocahzed basis. It is 
evident from Eq. ([5]) that if, for example, A = A, A2 = and the atomic mode S2 is on resonance with the field 
mode 1712- However, provided that {A,y4} > g, the modes si and nii will be off-resonant and will not contribute 
significantly to the dynamics on short time scales. Away from the resonances at A = ±A, both atomic modes are 
detuned from their corresponding field modes and the system dynamics becomes dispersive in nature. 

We consider three limiting cases: large hopping {A 3> {A, 5}), large detuning (A ^ {A,g}), and near-resonance 
{A + A ^ {A — A,g}). For each, we determine the effective dispersive-case Hamiltonian, which illustrates the 
important physical processes, and then calculate the explicit single-excitation dynamics from the full solution given 
in Eq. ®. 



III. LARGE HOPPING 



The first limit we examine is the large-hopping limit, in which A 3> {A,^}. The atoms are highly detuned from 
both delocahzed field modes ({Ai, A2} 3> <?). Transitions in which an excitation is passed from atom to field (or vice 
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versa) involve the exchange of a large amount of energy and so are strongly supressed. We can simplify the dynamics 
by constructing an effective Hamiltonian in which such transitions are absent. 
We begin by switching to the interaction picture with respect to 
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{bJa |e)j (e| +ujfa}-aj^ + A {a\a2 + a^di^ 



In this picture, the Hamiltonian becomes 



h.c. 



(11) 



(12) 



Since Ai and A2 are both taken to be large in this limit, the terms in this interaction Hamiltonian rotate rapidly in 
time. As in the standard rotating- wave approximation [3l|,[3l|, the time average of such rapidly rotating terms may 
be set to zero, thus eliminating the highly off-resonant transitions between atomic and field modes. 

To simplify the Hamiltonian in this manner, without entirely eliminating the atom-field interaction terms, we follow 
the method of James and Jerke |33|]. The authors define an effective Hamiltonian Tie// from the time derivative of 
the time-averaged evolution operator, 



where the time average of the operator 0{t) is defined as 

/oo 
f{t~t')6{t')dt'. 
-00 

Here, the probability density function f(t — t') acts to remove the high frequency components of the operator. 
By comparing Eq. p3p with the Schrodinger equation for the non-time-averaged evolution operator, 

i^U{t) = HilJ{t), 



(13) 



(14) 



(15) 



we can write the effective Hamiltonian as a series in Hi and U . For a system with a weak interaction, we discard 
terms beyond second order in Hj. The Hermitian part [37| of the effective Hamiltonian may thus be written as 



^e// - Hi{t) 



1 



Hj{t),Ui{t) - Hi{t),Ui{t) 



(16) 



The term Ui{t), given by 



Ui{t) = - / Hi{t')dt' 



(17) 



is the first order term in the expansion in Hj{t) of the evolution operator U(t). 
Applying this method to our system, we replace Ui{t) in Eq. (fT6|) by 



V{t) 



a; 



h.c. 



where iV{t) is the time integral of the interaction Hamiltonian. 

The effective Hamiltonian can thus be determined by taking the time-average of (|12p , , and 



Ai 



TO1S1+, rh\si 



A, 



171232+, rhls2- 



1 

a7 



1 

a; 



1 

a7 



1 

A~2 



m2S2+,m|si- 



(18) 



(19) 
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In the present limit, those terms rotating at rate Ai, A2, or 2A are considered to be fast, so their time averages may 
be set to zero. The effective Hamiltonian is thus given by the first two terms of Eq. divided by 2. Returning to 
the Schrodinger picture, the effective Hamiltonian becomes 



2 



■,2 



H' ^Hq+ 5^ (m,ml\e)j(e\ - mlm,\g) , {g\^ + GA {0^+02- + o^+o^-) , (20) 



where 



This Hamiltonian has three contributing terms. The first, _ffo, is the energy of the bare atomic and delocalized field 
modes, given by Eq. (|f f p . The summation represents a dispersive interaction between the atoms and the delocalized 
cavity modes m\ and m^. As with dispersive interactions between atoms and single cavities, there is no transfer of 
energy; the atoms, however, experience a Stark shift, the magnitude of which depends on the population of the field 
modes. Finally, the last term describes a transfer of energy between the two atoms without excitation of the field 
modes. 

The behavior of two atoms coupled to two highly detuned, non-interacting field modes described here closely 
resembles the behaviour of two atoms coupled to a single detuned field mode [ll| . In the single- mode case, the atoms 
are Stark-shifted by the field mode, with the shift proportional to g^/A (where A is the detuning between the atom 
and the field mode). In the two-mode case, the atoms experience a shift due to each mode, proportional to g^/A^. The 
total effect is the average of the two single-mode shifts. In both cases the interatomic energy transfer occurs without 
excitation of the cavity modes. The rate of transfer, however, differs between the two models. In the one-mode model, 
the rate is inversely proportional to the atom-field detuning. In the two-mode model, it is inversely proportional to 
the product A1A2 of the detunings between the two field modes and the atoms, and proportional to the intercavity 
hopping A. 

To find a more precise description of the system's behavior, we turn to the equations of motion given by Eq. ([5]). 
Taking the limit of large hopping, we obtain 

a(t) = e-'l'^i-GA)* j^(Q) (^^^ ^ -^(Q) sin(^t)] 
h(t) = e-^l-'z-GA)* [_ia(o) sin(yli) -I- 6(0) cos {At)\ 

C{t) = e-^l-'a + GA)* J^(Q) _ -^(Q) 

d{t) = e-^('^<>+GA)t [_ic(0) sin {GAt) + d{0) cos {GAt)] . (22) 

As expected, there is no transfer of energy from field to atom or vice versa. If the excitation is initially in one of the 
field modes, it will simply oscillate between the cavities at rate A; this is the same as the rate of intercavity transfer 
in the absence of atoms. If the excitation is initially in one of the atoms, it will oscillate between the two cavities 
at rate GA, without populating the field modes. Thus, at times t = {n + ^) n/GA, an excitation initially in atom 
1 will be fully transferred to atom 2. Increasing the atom-cavity detuning A increases this rate of transfer, as doing 
so decreases the product A1A2 of the detuning between the atoms and the two delocalized field modes, and thus 
increases G. Increasing the hopping strength A causes a quadratic decrease in G. This is sufficient to overcome the 
linear scaling of the transfer rate GA with A; the transfer rate therefore diminishes with increasing hopping strength. 

Figure [1] shows the evolution of the system in the large hopping limit, calculated numerically from the full Hamil- 
tonian in Eq. ([T]). The initial state has a single excitation localized in atom f . In agreement with the approximate 
solution of Eq. ([22]) . the probability of finding the excitation in atom 2 undergoes sinusoidal oscillation at a rate 
~ GA, and the probability of finding the excitation in either cavity mode is almost zero. 

Having seen how a general state in the single-excitation subspace behaves, we turn to the problem of transferring 
a qubit state from atom 1 to atom 2. The initial state is that given in Eq. ([S]), which consists of a superposition 
of zero- and one-excitation states. As the zero-excitation state is unaffected by the Hamiltonian, we need only 
consider the evolution of the single-excitation part. In terms of the general state of Eq. ([8]), the initial conditions 
are a(0) = 6(0) = d{0) — 0; the condition for successful state transfer at time r is (i(r) = c(0). This places two 
requirements on the interaction time t. First, to ensure that the final state of atom 2 has the correct ratio between 
the excited and ground states, we need GAt = (2?! -I- l/2)7r {n — 0, 1, 2, ... ). Second, r must be chosen such that 
{uja + GA)t = (2m — l/2)7r (m = 0, 1, 2, . . . ) in order for the phase difference between the two states to be correct. 
Both can be achieved if the system parameters are selected so that 

^a + GA _ Am -I 

GA ~ 4n -f I ■ ^ 
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FIG. 1: (Color online) Dynamics of atomic and field excitation in the case A, calculated by numerical solution 
of the full Hamiltonian, Eq. ((!]). The initial state has a single excitation in atom 1. The dashed red curve gives the 
probability, as a function of time, of finding the excitation in atom 2. The solid blue curve corresponds to the 
probability of excitation of the field modes. The system parameters are A = IQg, A = O.lg, and ujf ~ lOOOg. 



If only the first condition can be met, then, provided that the interaction time is accurately known, a single qubit 
operation can be applied to atom 2 to correct the phase difference and complete the state transfer. 

This transfer of states between atoms in separate cavities without populating the cavity modes has significant 
practical benefits. The method of transferring the atomic state to the field and then from the field to the second 
atom [3| is susceptible to decoherence from photon loss in the cavities or optical fibers. Qubit-qubit excitation transfer 
without populating the field modes, as in the proposal by Zheng and Guo [ll[ and the realization by Majer et al. p^ . 
minimizes decoherence due to photon loss. Furthermore, it is possible to extend this technique to decoherence- 
free entanglement generation, as suggested in [llj. For an initial state with a single excitation in atom 1 [a(0) = 
6(0) ~ c?(0) = 0, c(0) = 1], Eq. shows that the two atoms will be maximally entangled at times t such that 
T = (n + i) 7r/2GA. 



IV. LARGE DETUNING 



We next consider the limit of large detuning, in which A ^ {A,g\. As before, the atoms are highly detuned from 
both delocalized field modes, so atom-field transitions may be eliminated. The analysis proceeds similarly to that 
in the large hopping limit. This time, however, A is a small quantity and thus only those terms that rotate at rate 
Al or A2 are eliminated from the time averages of Eqs. (fT2|) . (fT8|) . and (fT9|) . The resulting Hamiltonian is, in the 
Schrodinger picture, 

H' = Ha+^GA {ajd]\e)j{e\ - a]aj\g) ^{g]^ + ^ (^rhlm2 - m\rhi) {^iz + ^2z) + GA{ai+a2- +CT2+<ti_) . (24) 

The atom-atom excitation transfer, given by the final term of Eq. (|24p . is identical to that found in the previous 
limit. The Stark shift terms, however, take a different form. In the large-hopping limit, the shift for each atom was 
determined purely by its interaction with the two delocalized field modes. In the current limit, with the atom-cavity 
detuning much greater than the hopping strength, the shift splits into two contributing factors. The first of these, 
proportional to GA, is due to the local interaction between each atom and the field mode contained within its cavity. 
(Note that ii A = the Hamiltonian reduces to that of two isolated far-off-resonance atom-cavity systems.) The 
second contribution, proportional to GA/2, is due to an interaction between the atoms and the delocalized field 
modes. The Stark shift due to this term (which, as ^ <C A, is much smaller than the local effect) is proportional to 
the difference between the populations of the two field modes. 
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The evolution of a single-excitation state reduces to the same equations of motion as in the limit A ^ A, given by 
Eq. (HH). As before, there is no transfer of energy between the field modes and the atoms, and an excitation that is 
initially in an atom will oscillate between the two atoms at rate GA. In contrast with the previous limit, increasing 
the hopping strength decreases the product A1A2, thus increasing both G and the interatomic excitation transfer 
rate. Conversely, increasing A decreases the transfer rate. 



V. COMPARABLE HOPPING AND DETUNING 



Finally, we turn to the limit in which the atoms are nearly resonant with one of the delocalized field modes. Unlike 
the previous two limits where the atoms were highly detuned from both modes, direct transfer of energy takes place 
between the atoms and the field, creating a different mechanism for interatomic energy transfer. Specifically, we 
consider the limit A + {|A — A|, g}, with A > 0. In this case the atoms are nearly resonant with the field mode 
7712 but highly detuned from mi. 

For simplicity, we do not consider the equivalent parameter regime for negative detuning, i.e. A — A 3> {|A + ^|,(7}. 
In this limit, the atoms are nearly resonant with mi and highly detuned from m2. Unlike our previous work on 
the ground state of the coupled cavity system, which showed a marked asymmetry between positive and negative 
detuning '23'|, the dynamics depends only on the absolute value of the detuning and thus the positive- and negative- 
detuning cases exhibit the same behavior. 

We derive the effective Hamiltonian by eliminating terms rotating at 2A and Ai from the time averages of Eqs. 
(fT2)) . (|18p . and (|19p . The near-resonance condition means that A2 is small and the corresponding terms must be 
retained. Returning to the Schrodinger picture, we find 



H' — Hq 
^2 



2Ai 



mjs2- + m2S2^ 

(25) 

mimj (|e)i(e| 4- |e)2(e|) - mimi (|.g)i(.g| + \g)2{g\) - (o-i+a'2- + cr2+(Ti- 



The second term in this Hamiltonian describes the transfer of energy between the near-resonant atomic and field 
modes S2 and m2. The third term represents the dispersive interaction between modes Si and mi. As in the previous 
limits, this is characterized by a Stark shift of the atoms (which depends on the population of mi) as well as a 
direct transfer of energy between the atoms without stimulating the field modes ^38j]. As g ^ g^/Ai, the evolution is 
dominated by the near-resonant term in Eq. ((25|) . so the dispersive term may be neglected for times t ^ Ai/g'^. 

To study the dynamics in detail, we turn to the equations of motion for the delocalized modes given in Eq. ([9]) . We 
apply the current limit, A + A'^ {\A — A\, g}, and make the additional assumption that modes S2 and m2 are very 
nearly on resonance, i.e. |A — yl| ^ g. If the field modes are initially unpopulated, i.e. ctj{0) = 0, the evolution of 
the atomic modes is given by 



f3i{t) ^ e-'^^-^+s /'^^^*/9i(0), 
/32(i) = e-^('^»-^^/2)t^^(Q)^Qg 



(26) 



The population of atomic mode S2 , represented by coefficient f32 (t) , oscillates at rate ~ g due to the excitation transfer 
between S2 and field mode m2. Mode si, represented by (3i{t), has a constant population because it is coupled only 
dispersively to the field. The phases of the two modes rotate at different rates, with the difference S between them 
given by 

The second term in this expression corresponds to the Stark shift of atomic mode si due to its dispersive interaction 
with field mode mi. 

Let us consider the case of exact resonance, A2 = 0. If the effect of the dispersive processes is negligible, i.e. 
g'^/Ai <C Tr/t, so that 5 ~ 0, the two delocalized atomic modes rotate at approximately the same rate. The behavior 
of the localized atom and cavity modes, in this case, reduces to simple sinusoidal oscillation. Taking the initial 
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FIG. 2: Dynamics of atomic and field excitation in the resonant case A = A, calculated by numerical solution of the 
full Hamiltonian, Eq. ([1]). The initial state has a single excitation in atom 1. The dashed red curve gives the 
probability, as a function of time, of finding the excitation in atom 2. The solid blue curve corresponds to the 
probability of excitation of the field modes. The system parameters are A — lOQg, A — lOOg, and uif = lOOOg. 

conditions a(0) = b{0) = d{Q) = 0, we find that 

a(i) = 6W = -^e— °*sin(gi), 
c(t) = c(0)e--"*cos2 (^^y 

rf(i) = 0(0)6--"* sin2(^^^, (28) 

in agreement with the results of Nohama and Roversi [sot . Excitation transfer between atoms 1 and 2 occurs at times 
T — (2n + l)TT/g with rt = 0, 1, 2, . . . . In contrast with previous limits, an excitation that is initially in an atomic 
mode will be passed to the field. This direct transfer of energy from atom 1 to the field, then from the field to atom 2, 
is responsible for the high rate of interatomic transfer. Though the transfer occurs much faster than in the dispersive 
limit, the speedup comes at the cost of exposing the system to decoherence via cavity loss. 

Figure [5] shows the numerically determined evolution of the system for the case of exact resonance. The probability 
of finding an excitation in atom 2 oscillates between approximately zero and unity at rate g. We note, however, that 
this probability does not return to exactly zero. This is due to the Stark shift of mode si. As discussed below, this 
shift becomes significant for large interaction times and causes beating in the evolution of the populations of atoms 1 
and 2. Additionally, in contrast with the previous limits, the field modes become populated; the probability of finding 
a photon in either cavity mode oscillates between zero and 0.5. 

Perfect state transfer requires that ^(t) = c(0), introducing an additional constraint on the interaction time r. As 
well as selecting r so that r = (2n + l)7r/g, it is also necessary to have cJqT — 2m7r. As shown in [30], these can be 
simultaneously achieved by setting uja = 2lg where I is an integer. If this is impractical, the interaction time may 
be chosen such that r = (2n + l)Tr/g, so the final state of atom 2 has the correct population; the phase difference 
between excited and ground states can then be corrected by means of a single-qubit operation. 

Away from exact resonance, S ^ and the phases of the two delocalized atomic modes no longer rotate at the 
same rate. The competition between the resultant oscillating phase difference and the excitation transfer between 
S2 and 7712 leads to beating behavior in the evolution of the local atomic modes. Again taking the initial conditions 
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FIG. 3: (Color online) Probability of excitation transfer in the near-resonant case, calculated by numerical solution 
of the full Hamiltonian, Eq. ([T]). The initial state has a single excitation in atom 1. The curve shows the probability 
of finding the excitation in atom 2. The system parameters are A = 100. Ig, A = lOOg, and Wf = lOOOg. 



a(0) = 5(0) — d{0) — 0, we obtain from Eq. ([9]) the expressions 

ic(0) 



a{t) ^ b{t) 



g-j(tJ+Ai/2)t 
C{t) = <, COS 



A 



dit) 



i(w+Ai/2)t 



+ 



A 



^2 



-iSt 



-iSt 



(29) 



Figure [3] shows the probability of transfer when the atoms are detuned from field mode m2 by O.lg. The fast oscillation 
occurs at a rate ~ g, as in the resonant case. For the first Rabi oscillation, the probability of transferring an excitation 
from atom 1 to atom 2 is nearly unity. The effect of the detuning is to reduce the transfer probability over subsequent 
cycles; the maximum probability oscillates at rate d. 

This beating can be avoided for certain values of the system parameters. The difference 6 between the rates of 
rotation of the atomic modes si and S2, has two contributing factors: the detuning between S2 and m2, and the Stark 
shift of mode si. If these two factors have equal and opposite effects so that S = 0, the evolution of the system will 
reduce to the simple oscillation described by Eq. [28l Conversely, beating will occur even in the resonant case on 
timescales that are sufficiently long for the dispersive processes to become significant, i.e. t ~ l/<5 = Ai/g^. 



VI. DISCUSSION 



We have studied the dynamics of a system consisting of two cavities, each interacting with a two-level atom with 
coupling strength g and detuning A. The cavities are coupled through a photon hopping term with strength A. The 
electromagnetic field in such a system may be described either in terms of local cavity fields or delocalized modes. 
The behavior of the system is governed by the relative values of the three parameters g, A, and A. Let us assume that 
the values of g and A have already been set. Three different regimes of behavior may then be identified, depending 
on the value of A. When ^ 0, the two atom-cavity systems evolve independently. In the limits A <C A and A^ A, 
the field may be eliminated from the dynamics and the atoms interact directly. When A = ±A, the atoms interact 
resonantly with one of the delocalized field modes. 

This variety of interaction, governed by the single parameter A, could allow coupled cavities to form a useful and 
highly adaptable element of a QIP system. By selecting different values for the atom-cavity detuning and the hopping 
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strength, it is possible to choose the properties of the atomic state transfer. In the dispersive hmit, the transfer will 
take place slowly, but since the field modes are never excited the process is protected against decoherence. When 
the atoms are on resonance with one delocalized field mode, the transfer is more rapid, although it is subject to 
decoherence due to cavity loss. 

Such dispersive and resonant methods of state transfer could also be realized in a system of two atoms interacting 
with a single cavity. Coupled cavity systems, however, offer additional advantages. Qubits coupled to separate cavities 
could be more easily isolated from each other, eliminating cross-talk. If the cavity-cavity coupling is turned off, the 
two qubit-field systems become independent, and thus able to act as part of separate QIP subsystems. Furthermore, 
by controlling the hopping strength, the interaction between the qubits can be controlled without manipulating the 
qubits themselves. 

There are several possibilities for experimental realizations of coupled-cavity systems. Within the cavity QED 
setting, high Q-factor cavities could potentially be coupled by means of an imperfectly reflecting mirror. Atoms 
would be either trapped in the cavities or passed simultaneously through them. Another potential setting is that 
of semiconductor microcavites. A system of two coupled microdisk resonators, each with an embedded quantum 
dot, has already been demonstrated 34]. Perhaps the most promisingsystem, though, is circuit QED. Experimental 
work on single qubit-stripline systems is already well established [3, Hj. Majer et al. have demonstrated that it is 
possible to couple two superconducting qubits via a single resonator without significantly populating the resonator 
modes. Furthermore, a scheme for controllable coupling between stripline resonators has been proposed recently [ssj . 
A pair of qubit-stripline systems, connected by a controllable coupling element, would provide a flexible and realistic 
experimental implementation of the coupled-cavity model. 
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